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Do phenomenological master equations with memory kernel always describe a non-Markovian 
quantum dynamics characterized by reverse flow of information? Is the integration over the past 
states of the system an unmistakable signature of non-Markovianity? We show by a counterexam- 
ple that this is not always the case. We consider two commonly used phenomenological integro- 
diflferential master equations describing the dynamics of a spin 1/2 in a thermal bath. By using a 
recently introduced measure to quantify non-Markovianity [H.-P. Breuer, E.-M. Laine, and J. Piilo, 
Phys. Rev. Lett. 103, 210401 (2009)] we demonstrate that as far as the equations retain their phys- 
ical sense, the key feature of non-Markovian behavior does not appear in the considered memory 
kernel master equations. Namely, there is no reverse flow of information from the environment to the 
open system. Therefore, the assumption that the integration over a memory kernel always leads to a 
non-Markovian dynamics turns out to be vulnerable to phenomenological approximations. Instead, 
the considered phenomenological equations are able to describe time-dependent and uni-directional 
information flow from the system to the reservoir associated to time-dependent Markovian processes. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc 



I. INTRODUCTION 

The study of non-Markovian open quantum systems 
has attracted extraordinary attention and efforts in re- 
cent years Many analytical methods and numerical 
techniques have been developed to treat non-Markovian 
processes In addition to their importance in ad- 

dressing fundamental questions [12,] . this is mainly due 
to the applications non-Markovian systems find in many 
branches of physics. Non-Markovian processes appear in 
quantum optics 0, [l^, [ij] , solid state physics [l5[ , quan- 
tum chemistry [l6| , quantum information processi ne; jl7l| , 
and even in the description of biological systems [1^ . 

Recently, several more rigorous definitions and quan- 
tifications of non-Markovian behavior in open quantum 
systems have been proposed [l^-[2^. In fact, in the past 
the concept of non-Markovian dynamics has been quite 
loosely defined. The term non-Markovian process could, 
e.g., stand for: Not describable by a master equation with 
Lindblad structure, or leading to non-exponential decay, 
or characterized by a time-dependent generator, or in- 
volving an integral over the past states of the system. 

Quite often it has been argued that the treatment of 
non-Markovian dynamics necessarily requires solving an 
intcgro-differential equation for the reduced density ma- 
trix of the system. However, it has been shown that 
master equations which are local in time can also repre- 
sent the memory effects of a non-Markovian process (see, 
e.g., Ref. [H and references therein), without the need to 
take into account a time integration over the past his- 
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tory of the system. Here we go a step forward: Not 
only that memory kernel master equations are not the 
unique tool to treat non-Markovian dynamics, but we 
also demonstrate that the presence of a memory kernel 
alone does not guarantee the non-Markovian character 
of the process associated to the reverse flow of informa- 
tion from the system to the environment. This surprising 
result is obtained by applyinga recently proposed mea- 
sure of non-Markovianity [20l, ',2111 to two quite commonly 
used non-local master equations: The generalized mem- 
ory kernel master equation discussed by one of us and 
the post-Markovian Shabani-Lidar master equation flo|, 
both used to study the time-evolution of a spin 1/2 in 
a thermal bath. Our results are connected to the issue 
of phenomenological vs. microscopically derived master 
equations in quantum optics which do not always pro- 
duce coinciding results in all of the relevant parameter 
regimes [IJj. Here we show that there can be also qual- 
itative differences in addition to the quantitative ones 
when the two approaches are used. 

Before proceeding with our treatment we would like to 
emphasize that we do not intend to discourage from the 
use of memory kernel master equations, which indeed in 
many cases constitute a fundamental tool to study non- 
Markovian systems (see for example Ref. Instead, we 
rather want to point out that the mathematical features 
of phenomenological kernel equations do not guarantee 
that the feedback effects from the environment and the 
physical characteristics of non-Markovian dynamics are 
taken into account. 

We begin by briefly describing the measure for non- 
Markovianity introduced by three of us in Refs. (20l[2l|. 
Then we present the master equations under investiga- 
tion and their solutions for the density matrix of the re- 
duced system. Finally, we use these solutions to evaluate 
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the degree of non-Markovianity for the two quantum pro- 
cesses, concludmg with a discussion and some remarks. 



II. MEASURE FOR NON-MARKOVIANITY 

The construction of the measure for the degree of non- 
Markovianity in open systems is based on the definition 
of Markovian processes as those that continuously reduce 
the distinguishabihty of quantum states [20l|. One can 
interpret this loss of distinguishabihty as a flow of in- 
formation from the open system to its environment. By 
contrast, in a non-Markovian process there exists a pair 
of states the distinguishabihty of which grows for certain 
times. This growth of the distinguishabihty of states can 
be interpreted as a reverse flow of information from the 
environment to the open system which is defined to be 
the essential feature of non-Markovianity [20', ^27\ . 

An appropriate measure for the distinguishabihty be- 
tween two quantum states given by density matrices pi 
and p2 is the trace distance j25| 



^Tr|pi -P2I, 



(1) 



where \A\ = V A'^A. The trace distance represents a met- 
ric on the space of physical states. It has the important 
property that all quantum operations, i.e., all completely 
positive and trace preserving (CPT) maps are contrac- 
tions for this metric. Given a pair of initial states ^1,2(0) 
the rate of change of the trace distance under the time 
evolution is defined by 



presence of the time-integral over the past history, the 
two master equations do not describe any feedback of 
information from the environment to the system and are 
thus memoryless in this sense. However, it is important 
to note that the treated master equations can describe 
time-dependent uni-directional flow of information from 
the system to the environment. 



III. MEMORY KERNEL AND 
POST-MARKOVIAN MASTER EQUATIONS 

We first present a paradigmatic example of a phe- 
nomenological memory kernel master equation, describ- 
ing the dynamics of a spin 1/2 interacting with a bosonic 
reservoir at temperature T under the rotating-wave ap- 
proximation. 



dpjt) 
dt 



k{t')Cp{t ~ t')dt' . 



(4) 



Here, p{t) represents the reduced density matrix of the 
spin, k[t) is a memory kernel function containing infor- 
mation about the properties of the reservoir, and C is 
a Markovian superoperator. This superoperator is given 

by 



Cp = (2cr_/3Cr+ - cr+cr_p- p(7+cr_) 

+ -^(2cr+pcr_ - o^a+p - p(T-(T+), 



(5) 



a(t,pi,2(0)) = -D(pi(0,P2W)- 



(2) 



A given process is said to be Markovian if for all pairs 
of initial states the rate of change of the trace distance 
is smaller than zero for all times, i.e., a{t, pi^2{t)) l£ 0. 
Thus, a process is defined to be non-Markovian if there 
exists a pair of initial states /3i,2(0) and a certain time 
t at which the trace distance increases, a(t, pi,2(t)) > 
0. As shown in Ref. [l^, [2l| one can construct on the 
basis of this definition a measure for non-Markovianity 
which represents a functional A/'($) of the corresponding 
quantum dynamical map $. This measure is defined as 
the maximum over all pairs of initial states of the total 
increase of the distinguishabihty during the whole time- 
evolution: 



7V($) = max 

Pl,2(0) 



a{t, p,,2{0))dt. 



(3) 



cr>0 



In the following we are going to prove that for the 
dynamics of a simplified spin-boson model generated by 
the generalized memory kernel master equation Q and by 
the Shabani-Lidar post-Markovian master equation [lO| 
the rate of change of the distinguishabihty of any pair 
of states is always negative, implying that the measure 
of non-Markovianity is equal to zero. Thus, despite the 



with 7o being the phenomenological dissipation constant, 
N the mean number of excitations of the reservoir, and 
a± the usual raising and lowering operators of the spin. 
We consider a widely used form for the memory kernel 
function, namely an exponential function 



k{t) = 76^'^*. 



(6) 



The memory kernel master equation (jJl can be solved 
using the method of the damping basis . In Ref. Q 
the solution for the components of the Bloch vectors was 
derived from which one easily obtains the solution for the 
spin density matrix p{t) corresponding to a generic initial 
state. Using the basis {|0), |1)} of the eigenstates of az 
the elements of p{t) can be written as 



pii{t) = uit)pn{0) + v{t)pQoiO), 

poo{t) = (1 - w(<))pii(0) + (1 - «(t))poo(0), 

pio(t) = z(t)/9io(0). 



(7) 



where u{t), v{t), and z{t) depend on the damping matrix 

A = diag{Ai, A2, A3} and the translation vector 7^ — 
{Ti,T2,T3) asM(i) = (l+r3+A3)/2, vit) = (l-(-T3-A3)/2 
and z{t) = Ai. The damping matrix elements and the 
translation vector components can in turn be expressed 
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as 







\l= \2 = ^MiR/2,t), 

Ti=T2 = 0, 

^ 2iV + 1 ' 
where the function S,M{R-i) is given by 



(8) 



-p-7t/2 



1 



: sinh 



cosh 



(9) 



withi? = 7o(2iV + l)/7. 

An interesting master equation which interpolates be- 
tween the generahzed measurement interpretation of the 
Kraus operators and the continuous measurement inter- 
pretation of the Markovian dynamics is the Shabani- 
Lidar post-Markovian master equation. The general form 
of this master equation is [loj 



dt 



k{t')exp{Ct')p{t-t')dt', 



(10) 



where once more p{t) is the density matrix of the reduced 
system, k{t) the Shabani-Lidar memory kernel, and C is 
the Markovian superoperator. In the following C is taken 
to be of the form of Eq. ([5]), and the kernel function 
is again an exponential function given by Eq. ([5]). The 
solution for the density matrix elements of the spin can 
be written again in the general form of Eq. ([T]), where 
u{t), v(t), and z{t) depend in the same way as before on 
the damping matrix elements and the translation vector 
components They in turn have the same analytic 
expressions of Eqs. ((8|) with the exception that £,M{R,t) 
has to be replaced by the quantity £,p{R, t) which is given 

by 



-/? H- 1 

£.piR,t) =cxp ( —-ft 



1 



V^-riR) 
cosh 



: sinh 



(11) 



with r{R) = AR/{R+lf and R = 7o(2Ar + l)/7. 

In Ref. i7;] the conditions for the positivity and com- 
plete positivity of the dynamical maps associated to the 
master equations ^ and (fTO|) were studied. There it 
was found that 4i? < 1 is a necessary and sufRcient con- 
dition for the positivity of the dynamical map associated 
to the memory kernel master equation, while complete 
positivity is satisfied only for moderate and high temper- 
atures of the reservoir. On the other hand, the dynamical 



map corresponding to the post-Markovian master equa- 
tion ([TU]) is always completely positive. These results 
are in agreement with Ref. [1| where it was noticed that 
the memory kernel master equation can be derived from 
the post-Markovian one in the limit in which the phe- 
nomenological dissipation constant 70 is much smaller 
than the reservoir correlation decay rate 7, suggesting 
that the post-Markovian master equation is somehow 
more fundamental than the former one. Therefore, while 
for the post-Markovian Shabani-Lidar master equation 
we can freely investigate the non-Markovianity for the 
whole range of parameters, in the case of the memory 
kernel master equation we are restricted to the condi- 
tions 4i? < 1 and > 1. 

Having constructed the solution of the two master 
equations we can now determine the rate of change of 
the distinguishability given by Eq. which leads to 



a{t)±a{t) + \b{t)\j-^m\ 



(12) 



where a{t) = p\^{t) - pl^{t) and b{t) = p\^{t) - pj^it) 
represent the differences of the populations and the co- 
herences of the density matrices pi (t) and p2 (t) . It is easy 
to see that these function are equal to a{t) = A3(pJ]^(0) — 
p?i(0)) = hao and b{t) = Ai(p}o(0) - pUO)) = A160, 
with A3 = ^M(p){R,t) and Ai = ^a/(p) (-R/2, i) for the 
memory kernel and the post-Markovian master equation, 
respectively. The derivative of the trace distance can thus 
be written as 



ait) 



^aiaR.tf + wm^ty 



(13) 

with f (i?, t) = ^M(p) {Ri t) in the two cases. At this point 
we need to study the properties of these two functions. 

Let us consider first the memory kernel master equa- 
tion: Under the condition 4i? < 1, ^M{R,t) is a posi- 
tive, monotonically decreasing function which means that 
f^iuiR.t) < and ^6/(^/2,0 < 0. Since and |5oP 
are obviously positive for any pairs of states we thus 
have a{t) < 0. We conclude that the rate of change 
of the distinguishability of any pair of initial states al- 
ways decreases and, hence, the flow of information from 
the system to the environment is never inverted during 
the dynamics and non-Markovian effects do not appear. 
Analogously, also in the case of the post-Markovian mas- 
ter equation ^p{R, t) is a positive, monotonically decreas- 
ing function such that (j{t) < 0. This implies that also 
the post-Markovian master equation does not describe 
any feedback of information from the environment to the 
open system. In contrast to purely Markovian evolution 
fulfilling the semigroup property, the dynamics generated 
by the memory-kernel and the post-Markovian master 
equations can be classified as time-dependent Markovian 
processes since they do not fulfill the semigroup property 
while the information flow is nevertheless uni-directional 
from the system to the reservoir. This is demonstrated 
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by the time-dependence of the decay rates in the corre- 
sponding time-local description given below. 

Indeed, both master equations studied here can be 
written in the time-convolutionless form 



dp{t) ^ 7i(t) 
dt 2 

, 72 (i) 



2 

73(0 



(2o'+p(T_ — a^a^p — pa^cr^) 
{2azpaz - a^a^p - pa^a^), (14) 



where 



72 (i) 



'2N+1 ^{R,t) ' 
2N+1 ^{R,t) ' 



(15) 



and ^{R,t) — ^M(p)(^ii) in the two cases. Thus we 
see that the integro-difFerential master equations given 
by Eqs. (H])-® and by Eq. (flUl) can be transformed into 
a form which is local in time and does not involve any 
time-integration over a memory kernel. 

The decay rate 73 (i) in Eq. (flSj) is always negative. 
It follows that the dynamical map $ corresponding to 



the master equation is nondivisible |21|. On the other 



hand, we have found above that the process is Markovian. 
Thus we have an explicit example of a nondivisible quan- 
tum process with zero measure for non-Markovianity, 
A/'($) — 0. The existence of such processes was already 
conjectured in Ref. [2T'|. Physically this means that the 
influence of the decay channel with a negative rate is over- 
compensated by the effect of the other channels with pos- 
itive rates, such that the distinguishability of quantum 
states is still monotonically decreasing. We include these 
non-divisible processes which have uni-directional infor- 
mation flow into the class of time-dependent Markovian 
processes. This class also includes processes whose de- 
cay rates are time-dependent and positive quantities [2l| . 
We emphasize that the time-dependent uni-directional 
(time-dependent Markovian) processes and the reversed 
information flow (non-Markovian) processes have impor- 
tant fundamental differences, as described recently, e.g., 
in Refs. [l[2l|2l|. 

Going back to the memory kernel master equation Q , 
it is also interesting to notice that the non-appearance of 
memory effects depends on the restrictions of the range 
of parameters imposed by the requirement of positiv- 
ity. In fact, when positivity breaks down for 4i? > 1, 
the hyperbolic sine and cosine of Eq. ^ are replaced 
by trigonometric sine and cosine. The function ^]\j(R,t) 
then shows damped oscillations, its derivative has no def- 
inite sign, and, consequently, there can be intervals of 
time in which the rate of change of the trace distance 



a{t) becomes positive implying that non-Markovian ef- 
fects appear. We mention that a violation of the posi- 
tivity of the dynamical map in phenomenological master 
equations was previously studied by Barnett and Stcn- 
holm in Ref. [27j]. There it was shown that the intro- 
duction of an exponential memory kernel function in the 
dynamics of a damped harmonic oscillator can lead to 
blatantly non-physical behavior. 



IV. DISCUSSION AND CONCLUSIONS 

We have applied a recently developed measure for 
the degree of non-Markovianity of quantum processes to 
the dynamical solutions of a simplified spin-boson model 
given by two widely used integro-difFerential master equa- 
tions. It has been demonstrated that, as long as the re- 
quirement of the positivity of the associated dynamical 
maps is fulfilled, no non-Markovian behavior occurs, i.e., 
the measure of non-Markovianity is equal to zero. This 
means that the phenomenological memory kernel mas- 
ter equations considered here are not able to describe a 
genuine non-Markovian behavior involving a backflow of 
information from the environment to the open system. 

Recently, the exact memory kernel master equation for 
a two-state system coupled to a zero temperature reser- 
voir has been constructed [2^, showing that in this case 
the structure of the master equation given by Eqs. (HI and 
([S]) is incompatible with a non-perturbative treatment of 
the underlying microscopic system-reservoir model. The 
perturbation expansion of the exact memory kernel re- 
veals that in higher orders a new decay channel appears 
in the superoperator ([5]) which is not present in the stan- 
dard Born approximation. Thus, while the exact mem- 
ory kernel master equation describes correctly all non- 
Markovian features of the model, approximation schemes 
and phenomenological models can lead to strong restric- 
tions in the treatment of non-Markovianity. 

Generally, one might be tempted to think that the 
introduction of a memory kernel necessarily leads to a 
dynamics with non-Markovianity and memory effects. 
However, our results demonstrate that one needs to be 
cautious when characterizing the physical properties of 
open systems only through the mathematical structure 
of their equations of motion. The presence of an integral 
over the past history in a phenomenological or approx- 
imate master equation does not necessarily guarantee a 
proper description of memory effects, namely, the feed- 
back of information from the environment to the open 
system. Even though memory kernel master equations 
certainly provide a very useful tool for the description 
of non-Markovian quantum processes, our results lead to 
the following questions: Which of the commonly used 
phenomenological or approximate memory kernel master 
equations are able to reproduce the key features of non- 
Markovianity? Can one formulate general conditions for 
the structure of the memory kernels which guarantee the 
presence of these features in the dynamics? 
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